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Abstract 

The general structure of the 1/m corrections at zero recoil is studied. The relevant 
matrix elements are forward matrix elements of local higher dimensional operators 
and their time ordered products with higher order terms from the Lagrangian. These 
matrix elements may be classified in a simple way and the analysis at the non recoil 
point for the form factor of heavy quark currents simplifies drastically. The second 
order recoil corrections to the form factor \ia\ of the axial vector current, relevant for 
the |V c b| determination from B — > D* decays, are estimated to be —5% < Hai — 1 < 0. 
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1 Introduction 



Heavy quark effective theory is by now the standard description for systems with one heavy 
quark The additional symmetries appearing in the limit of infinite heavy quark 

mass yield model independent relations between form factors appearing in the description 
of heavy hadron exclusive weak decays. Aside from that, heavy quark symmetries also 
yield statements about the normalization of some form factors at zero recoil, i.e. the point 
where the velocities of the initial and final hadrons are equal. This fact has a very large 
phenomenological impact; it allows e.g. to perform a model independent determination of 
\V c b\ by extrapolating to the endpoint of the lepton spectrum in the decay B — > D*iv. 

QCD radiative corrections as well as recoil corrections have been studied already to next- 
to- leading order [|Hj. While QCD radiative corrections may be studied systematically, the 



recoil corrections in general need new, non-perturbative input, which may be supplied, for 
instance, by model estimates. For the case of the determination of \V c b\ from B — > D*lv the 
leading recoil corrections vanish at the non-recoil point due to Lukes theorem |ll|] and the 



next-to-leading ones have been considered by Falk and Neubert [|T^], who parameterized the 
form factors to order 1/mg, also off the point of equal velocities. 

However, as will be discussed below, that the analysis of the 1 /mg corrections at zero 
recoil simplifies enormously, since then only forward matrix elements (i.e. matrix elements of 
operators between mesons moving with the same velocity) appear. In addition, the algebra 
of Dirac matrices simplifies and one may obtain a simple expression for the next-to-leading 
recoil corrections at the point of equal velocity. This expression involves forward matrix 
elements of operators of higher dimension and also time ordered products with higher order 
recoil terms from the Lagrangian. The expressions we obtain have a simple interpretation, 
but its numerical evaluation needs input beyond heavy quark effective theory. 

Recently, the methods of the heavy mass expansion have been applied also to inclusive 
decays by combining the method of operator product expansion with heavy quark effective 
theory |13[-||17||. This approach yields the heavy mass expansion for decay rates and also for 



decay distributions; the leading term in this expansion is the free quark decay rate and the 
corrections may be studied systematically. Of course, the higher order corrections need non- 
perturbative input, which is again forward matrix elements of higher dimensional operators 
and time ordered products of such operators with higher order terms from the Lagrangian. 
Thus the same matrix elements appear as in the model independent determination of V c b- 

Finally, the relation of the heavy hadron mass to the mass of the heavy quark is also given 
in terms of a 1/mg expansion. Higher orders are again given by forward matrix elements 
of higher dimensional operators needed as non-perturbative input to relate the heavy quark 
mass with the heavy hadron mass. 
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In the present paper, a systematic study is performed for these forward matrix elements 
appearing in all higher order calculations at zero recoil, including the relevant time ordered 
products with higher order terms of the Lagrangian. It turns out that all the forward matrix 
elements may be classified very simply and the relevant matrix elements for calculations up 
to order l/m% are given explicitly. 

The classification performed here allows to simplify the analysis of the recoil corrections 
to heavy quark weak decay form factors at v — v' enormously, compared to the case off the 
equal velocity point. As an application the analysis for weak decay form factors is performed 
at the non-recoil point up to second order in the heavy mass expansion for the case of b — > c 
transitions and our results are compared with the ones obtained by Falk and Neubert fT2|l . 

In the next section, a general discussion of the parametrization of the generic forward 
matrix element is given. It is split into three subsections. First we consider local higher 
dimensional operators in some detail and give numerical estimates for the matrix elements 
of operators up to dimension seven. In the second subsection we shall consider time ordered 
products with higher order terms from the Lagrangian. Finally, in the third subsection we 
consider the relation between the mass of the heavy meson and the mass of the heavy quark 
as a toy example, where the forward matrix elements play a role. 

The general discussion of the forward matrix elements is then applied in section 3 to the 

1 / vt?q corrections of the normalization of the weak decay form factors in the decays B — > Div 
and B — > D*£u. The relevant form factors h + for B — > Div and Ua\ for B — > D*lv decays 
are discussed at the non-recoil point up to second order in the heavy mass expansion. We 
include also a semi- quantitative analysis ad estimate the size of the corrections relevant for 
the V c b determination to order 1/rriQ. 

2 Higher Dimensional Operators and their Forward 
Matrix Elements 

Higher order terms in the heavy mass expansion of weak transition matrix elements originate 
in general from two sources. The first source is the heavy mass expansion of the operators 
for a heavy quark Q appearing in the weak transition Hamiltonian. At the matching scale 
ttiq, this amounts to the replacement 



Q(x) = e~ im ^ vx) 1 + J2 



( 



2m Q 



ivD 



k 




Q v (x) 



(1) 
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where v is the velocity of the heavy hadron and Q v is the operator of a static heavy quark 
moving with velocity v. Furthermore, it is convenient to define 

D^ = {9 liV -v^ v )D" vD ± = 0. (2) 

These terms in general lead to contributions, which are matrix elements of local operators. 

Secondly, also the Lagrangian of full QCD is expanded in I/tuq and the higher orders 
in I/rag are treated as perturbations. This leads to time ordered products involving these 
higher order terms of the Lagrangian and the weak transition operator. The corrections of 
higher orders in the 1 / tuq expansion to the Lagrangian are given at tree level by |§ 

oo oo / 1 \ 3 

Cj = E4 J) = Ek7 Q v (-i^)(ivDy-\i^)Q v (3) 
j=i j=i \ Q J 

In every order j it is convenient to split C^p into a generalized kinetic energy operator 
and chromomagnetic moment operator defined as 

= Q v {iDi){-ivDy-\iD La )Q v (4) 
gU) = -iQ v (iDi){-ivDy- l {iD^)a aP Q v (5) 

The interpretation of these time ordered product terms is obvious. The heavy hadron 
states of full QCD still depend on the heavy quark mass, and this dependence is also treated 
in a 1/rriQ expansion. The leading term is the state taken in the infinite mass limit, the 
"static state", which is the convenient one for a heavy quark effective theory calculation, 
since it does not depend on the heavy mass any more. The time ordered products account 
for correct mass dependence of the full QCD state, and the matrix elements then have to be 
evaluated using the "static", mass independent state. 

2.1 Local Operators of Higher Dimension 

The generic operator of dimension n + 3 appearing in the contexts mentioned above is of the 
form 

°21«-*. = QviiD^iD^) ■ ■ ■ {iD, n )TQ v (6) 
where V is an arbitrary Dirac matrix. 

The Dirac matrix appearing in (|6]) may be expanded into the 16 basis Dirac matrices 1, 
75) Iij, 757^ an d cr^v However, the matrix T is sandwiched between projectors 
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which are contained in the heavy quark fields Q v . This projection amounts to the replace- 
ments 

1 — P + = \{l + i>) 7,. — = v»P+ (7) 

7^75 — ► P + 1^>P+ = s M 7s — ► P+I5P+ = (8) 
(-i)<V — ► P+{-i)(T l)V P + = iv a e ailv ps (9) 

where we have defined the spin matrices s M , which are the generalizations of the Pauli 
matrices for the frame moving with velocity v. They satisfy the relations 

s^s u = (-g^ u + v^v u )P + + ie afMU(3 v a s (S vs = (10) 

Consequently, the Dirac matrix T sandwiched between the projectors may be expanded into 
the four matrices 1 and 

P+TP+ = l -P + Tr {P + T} - l -s, Tr , (11) 

and it is sufficient to consider only the two operators 

°tL-» n = QviiD^iiD^) ■ ■ ■ {iD, n )Q v (12) 
^:„,.„.,a = QviiD^iiD,*) ■ ■ ■ (iD^s.Qv (13) 

In the following we shall consider the matrix elements between the ground state pseu- 
doscalar and vector mesons. There are two different cases to be studied. In the first case the 
initial and the final state are both either 0~ or 1~; in the second case the initial state is 0~ 
and the final state is 1~ or vice versa. All these different cases are related by heavy quark 
spin symmetry, which implies the relations 

Q v \H(v))=>y 5 fiQ v \H*(v,e)) Q v \H*(v, e)> = l5j tQ v \H(v)). (14) 

where \H(v)) and \H*(v,e)) denotes the 0~ and \~ ground state meson respectively. 

Further restrictions on the structure of the matrix elements may be obtained from the 
equations of motion for the heavy quark and the gluons 



ivDQ v = [(iDj, (iD u )]] ab = Arra s £ 

where the sum runs over all quark flavors and a, b, d are color indices 



1 A _ 



(15) 
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We shall start the discussion with the case where both, initial and final state are either 
0~ or 1~. Spin symmetry relates the matrix elements of 0~ mesons with the ones of the 1~ 
case in the following way 

(H(v)\0^ 2 ...jH(v)) = (H*(v,e)\OV^...jH*(v,e)) (16) 

mv)\0^...^\H(v)) = -l(H%v,e)\Ol% 2 ...^H*(v,e)), (17) 

and we shall consider in the following only the 0~ case. 

The equations of motions for the heavy quark imply that the matrix elements of both 
and have to vanish, if the first or the last index, i.e. /ii or /i n is contracted with 
the velocity v. Furthermore, the spin vector s is also orthogonal to the velocity and thus the 
matrix elements have to satisfy 

v^(H(v)\0^ 2 ...JH(v)) = v^(H(v)\0£1 2 ...JH(v)) = 
v^mv)\0^„ n , x \H{v)) = v^(H(v)\Ol s l n , x \H(v)) = 
v\H(v)\0^l 2 ...^ x \H(v)) = (18) 

Note that contractions with any other index may be related to a gluon field strength 
[(iDp), {iD u )\ = igG^, e.g. 

{ivD^iD^Q^Hiv)) = -igv a G aiXn Q v \H{v)), (19) 

and are thus in general nonzero. 

Combining the information form the spin structure and the restrictions form the equation 
of motion of the heavy quark one obtains for the forward matrix element of the general 
expression 

(H(v)\Q v (iD a )(iD vl ) ■ ■ ■ {iD Vn _ 2 ){iD p )Q v \H{v)) = 2M H [g a(3 - v a v p ]A vl ... Vn _ 2 (20) 

The tensor A is constructed from g^. and v^. It is a simple combinatorical exercise to 
show that the number N of independent scalar parameters is 

N/2-1 n k 1 

" H = 1+| -" S (i) HCT (21) 

where n > 2 and [n] = n for n even and [n] — n — 1 for n odd. The number of independent 
parameters grows rapidly, the first few are N(A) = 2, N(5) = 4, N(6) = 13, N(7) = 41 and 
N(8) = 196. 
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The matrix elements of 0^ s > are parity odd quantities. The general form of these matrix 
elements, which is compatible with the restrictions (|T8|), is given by 



[H(v)\O ttltl ... lin _ 2 p.x\H(v)) = 2M H dme uapx v u B Ul ... Un _ 2 (22) 



+ 2M H d H \g a p-v a vp)C^ 



Vi — "n-2;A 
(2) 



+ 2M H d H [g a \-v a v x ]C ( i ;L Un 
+ 2M H d H [gpx-vpv x ]C^ 



P 

v\—v n -i,et 



where dn = 3 for a pseudoscalar meson and cf# = — 1 for a vector meson. The tensors C^ k > 
are parity odd and vanish, if the last index in contracted with v. 

Up to dimension seven the number of parameters is still manageble, and some of them 
are more or less well known numerically. The only nonvanishing matrix element between 
heavy meson states of the dimension 3 operators is 

(H(v)\Q v Q v \H(v)) = 2M H (23) 

and its value is given by the choice of the normalization. Here M# is the mass of the heavy 
meson in the static limit. 

All matrix elements of the dimension 4 operators vanish due to the equations of motion; 
all matrix elements of the dimension 5 operators are given in terms of two parameters Ai 
and A2 

(H{v)\Q v {iD a ){iD p )Q v \H{v)) = 2M H [g a(3 - v a v^\ x (24) 
(H(v)\Q v (iD a )(iDp)s x Q v \H(v)) = 2M H d H ie 1/a/3x v l/ ^\ 2 (25) 



where the prefactors are chosen to comply with the definition in [12]. The parameter A2 



corresponds to the leading term in 1/ rag for the mass splitting between the ground state 1 



and 0~ mesons ||18|| , while the kinetic energy parameter Ai is not related in an easy way to 
a measurable quantity. From QCD sum rule analyses one obtains values of Ai = —0.6 ±0.1 
GeV 2 fni , but these calculations have been criticized recently and a much lower value of Ai 



has been suggested using an improved sum rule technique |20| . On the other hand, bounds 
have been derived in a quantum mechanical framework indicating that Ai < —0.18 GeV 2 



2jJ . In the numerical studies presented below we shall vary Ai in some range and hence we 



shall use the values 

- 0.3 GeV 2 < Ai < -0.1 GeV 2 A 2 = 0.12 GeV 2 . (26) 
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The parameter A2 is scale dependent and we define A2 = A2(m&). 

The matrix elements of the dimension 6 operators are also given in terms of only two 
parameters p\ and p 2 

(Hiv^QviiDaXiD^iiD^KlHiv)) = 2M H [g a p - v a v p \v^ Pl (27) 

{H(v)\Q v (iD a )(iDp)(iD p )s x Q v \H(v)) = 2M H d H ie uaf3X v u v l ^p 2 . (28) 

In order to estimate pi we may employ the equations of motion for the gluon fields and relate 
this parameter to a forward matrix element of a four fermion operator 



- AM HPl = 4vra s 



(Qv, a Qv,b)(qbi>q a ) - -^-{Qv, a Qv,b){qbfqb 



\H(v)) (29) 



where a and b are color indices. Contracting ( [29D with v we may rewrite the matrix element 
as an interaction between two vector currents. Using the Fierz theorem we rearrange the 
quark fields in order to apply vaccum insertion, after which one is left with matrix elements 
of heavy light operators between the heavy meson and vacuum. These matrix elements are 
all related to the heavy meson decay constant fjj due to heavy quark spin symmetry. The 
estimate for the parameter p 1 reads under these assumptions 

1 N 2 - 1 

* = 2™'-N*-f* MB - (30) 



A similar estimate has been performed in [BTf . 



However, (|30D has the usual problem of a matrix element after factorization. The original 
matrix element defining p\ is expected to have a different behavior under renormalization 
group transformations as the result after factorization. In other words, one has to define at 
which scale factorization is performed. We shall factorize the matrix elements at the scale 
rrib and thus use the following set of parameters: a s = a s (rrib) = 0.2 and M# = 5.28 GeV. 
Varying the heavy meson decay constant between 150 and 200 MeV we obtain 

( Pl )V 3 = (300 - 450) MeV. (31) 

This number is of the same size as e.g. A^ 2 = 350 MeV. 

Finally, the forward matrix elements of the dimension 7 operators may be written 
in terms of two parameters 77 and r 

(H(v)\O l £ lltafj \H(v)) = 2M H ^[g a(3 - v a v f3 ]{g lillX2 r] 1 - v^v^n) (32) 
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while the general form of the dimension seven operator is more complicated 

(H{v)\0 { :l^ ]X \H{v)) = -IMuduie^v^g^B^ - v^B 2 ) (33) 

+2M H d H C {1) [g Ql3 - v a v p ] 
+2M H d H C {2) [g aX - v a v x }e p ^ 2l3 

One may again apply the equations of motion for the gluon field to relate this to a matrix 
element involving the light quark current. In this way one obtains a relation of the form 

2M H (Ari + r)(g a(3 -v a vp) (34) 

= -ATra s J2( H ( v )\ ((iD a Qv,a)Qv,b)(qbipq a ) - -rj-((iD a Q V: a)Qv,a)(qbiuqb) \H(v)) 

This may again be estimated by using the Fierz theorem and vacuum insertion. After 
factorization, using 

(H(v)\(iD a Q v )^ l5 q\0) = 3Af H M H [g afM - v a v^ A = M H — m Q , (35) 

one obtains 

_/V 2 — 1 - 

A V + r = 6na s ^^Af 2 H M H (36) 

With the same of parameters and under the same assumptions as above one obtains the 
estimate 

(Ari + r) 1/4 = (700 - 950) MeV, (37) 

where A has been varied between 400 and 600 MeV. 

From these estimates it seems that the heavy mass expansion works quite well, at least at 
the non-recoil point. All the parameters up to dimension seven behave like the appropriate 
power of some small scale A ~ 200 - 500 MeV which means that the expansion in powers of 
A/rriQ indeed has coefficients of order unity. 

The second case to be studied are matrix elements with a 0~ meson in the initial state 
and a 1~ in the final state, or vice versa. These matrix elements do not introduce any 
new parameters, since they are related to the ones considered above by heavy quark spin 
symmetry. However, one has to be a little more careful in this case, because one has to 
rotate the spin of only one of the heavy quarks. The forward matrix elements, which are 
considered here, involve only one velocity sector of heavy quark effective theory, and spin 
symmetry is a symmetry holding separately in each velocity sector. In order to rotate only 
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the initial or the final state heavy quark spin, one has to choose in a first step two different 
velocities for initial and final state, perform the spin rotation of one of the states using (|16|) 
in the corresponding velocity sector, and afterwards take the limit v' — > v. In this way one 
obtains the relations 

{H{v)\Q v {iD^) ■ ■ ■ (iD^)Q v \H{v)) = -(H(v)\Q v (tD^) ■ ■ ■ (iD^)(se)Q v \H*{v, e)) (38) 
{H{v)\Q v {iD^) ■ ■ ■ (tD, n )Q v \H*(v, e)> = -{H(v)\Q v (iD^) ■ ■ ■ (tD, n )(se)Q v \H(v)), (39) 

relating the matrix elements of between two 0~ or two 1~ states to the ones of 
between 0~ and a 1~ state, and vice versa. 



2.2 The Time Ordered Products with the Lagrangian 

The second type of matrix elements appearing in an analysis of higher order I/tuq corrections 
at zero recoil are time ordered products of the local operators discussed above and the terms 
appearing in the heavy mass expansion of the Lagrangian. 

We shall first consider the case of two different flavors q v and Q v . The simplest terms 
are the two point matrix elements 



-i 



-i) 



d 4 x (H q (v)\T [q v (iD^) ■ ■ ■ {iD, n )Q v K,^\x)\ \H Q (v)) 
d'x (H q (v)\T \q v {iD^) ■ ■ ■ (iD^)Q v g®(x)] \H Q (v)) 



(40) 
(41) 



where the operators without argument have to be taken at x — 0. JCq and Qq are the kinetic 
and chromomagnetic terms for the quark Q as defined above. 

The spin structure of the simplest two point matrix elements may be analyzed in the 
trace formalism 



-t 



-t 



d 4 x (H q (v)\T [q v rQ v K%\xj\ \H Q (v)) = -A Tr {m{v)TM{v)} (42) 
d 4 x(H q (v)\T[q v TQ v g^\x)} \H Q (v)) = \b Tr {(- i )a a/3 M(v)P + (-t)a a(3 M(v)} 



where T is a general Dirac matrix, which is a linear combination of 1 and s^, and M(v) are 
the usual representation matrices for the heavy ground state mesons 



M(v) 



{j} + 1)75 pseudoscalar meson 

— (J) + 1)^ vector meson, polarization e 



(43) 



and the normalization is chosen according to (E3 
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The matrix a a p in the expression for the chromomagnetic moment operator appears only 
between projection operators P + and it is convenient to switch to a representation using the 
Pauli matrices (|9[). In this representation one has 



gti) = iv^ x Q v (iD a )(ivDy(iDp)s x Q v 
and we write for the second equation of ([4^) 

d 4 x(H q (v)\T[q v TQ v g^(x)} \H Q (v)) = -B Tr { lxl5 M[v)Ts x M(v)} 



(44) 



.-%) 



The representation in terms of the Pauli matrices is very useful, as soon as more than 
one insertion of a chromomagnetic moment operator appears, since the spin structure of 
products of chromomagnetic operators correspond to products of the spin matrices s which 
may be reduced using the relation ([3~0|) . For example, the product of two chromomagnetic 
moment operator insertion may be written as 

d 4 xd 4 y(H g (v)\T[q v Q v g$\x)g%\y)] \H Q (v)) = -Tr {T^M~(v)s a s>>M(v)} 
= -Tr {T^(-g a , + v a vjM[v)M(v)} - Tr {T^ie pa ^H{v)s v H{v)} (45) 
where T parameterizes the light degrees of freedom 



and one obtains 



T aP = l -T^\v a v p - g af3 ) + ~T^e 



(46) 



d A x d 4 y(H g (v)\T [q v Q v g^(x)g^(y)\ \H Q (v)) = 2M H (T^ + d H T^) 

In this way one may easily identify the spin symmetry conserving and spin symmetry vio- 
lating contributions of such products. 

The equations of motion also imply restrictions on the matrix elements of time ordered 
products [16, 22 1 . In principle, one obtains the same relations as for the local terms, for 
example 

(H q (v)\T [q v (tvD)Q v IC^\x)} \H Q {v)) = 

However, there may be an ambiguity depending whether the derivative acts on the T symbol 
or not. If one also takes the derivative of the step functions comming from the T symbol, 
then one obtains a local contribution of the form 



(H q (v)\T \q v (tD,)Q v IC^(x)} \H Q (v)) ~ t5\x)v,(H q (v)\q v (tD^(tvDy (tD ±a )Q v \H Q (v)) . 
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which may in general be reabsorbed into a redefinition of the T product. However, in the 
applications discussed below only the perpendicular components of the derivatives defined 
in (0) enter the expressions as e.g. 

(H q (v)\T [q v (t^)Q v }C^(x)} \H Q (v)) = 0, 

and hence there will be no contribution from such terms. 

Finally, the flavor diagonal case may be discussed by inserting first the correction terms 
for the Lagrangian of the quark q 



and then consider the case q = Q. In this case one has insertions in both lines, the one 
corresponding to q and to Q. When the masses are equal, both insertions are parametrized 
by the same form factor. However, the spin structure is different; in particular, the insertion 
of the chromomagnetic moment operator yields a Pauli matrix s to the right of V for Q , 
while s occurs to the left of V for q. Thus one obtains for the examples studied above 



-i) 



d 4 x(H Q (v)\T 
d 4 x(Ho(v)\T 



'Q v YQ v K^\x)] \H Q {v)) = -2ATr {m(v)TM(v)} (47) 
QvTQ v g®(xj\ \H Q {v)) = -B Tr { lxl5 M~(v){r , s x }M(v)} (48) 

where {, } denotes the anticommutator of the two Dirac matrices. 



2.3 Simple Application: The Heavy Meson Mass 



The mass of a heavy hadron may be expanded in inverse powers of the heavy quark mass. 
The lowest order terms of this expansion have been considered and one may extend this 
analysis to higher orders using the above discussion of the forward matrix elements. 

The relation between the heavy meson mass m# and the mass of the heavy quark is given 

by 



m H = M H - 



1 



2M, 



(H(v)\T 



H 



£/(0) exp (^—i J d 4 x Ci(x] 



\H(v)) 



(49) 



where M# = niq + A is the mass of the heavy hadron in the limit mq — > oo 
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The 1/rriQ expansion of the hadron mass is obtained by inserting the expression ([3]) into 
the time ordered product. Up to order I/tuq one finds 



Tfijj = m Q + A 



2M H 



E 

3=1,2 



2m c 



(H(v)\K®(ti)\H{v)) + (H(v)\g^(0)\H(v)) 



1 



-(-i 



2M H 

3 \ 



2m c 



d 4 x (H(v)\T [(/C«(0) + gW(0j) (K^(x) + g {1 \x))] \H(v)) 



+0{l/m i Q 



(50) 



The matrix elements appearing here are exactly of the type considered above. To order 
1/mg there are the two parameters Ai and A2, while to order 1/rrin one has not only local 
operators, but also time ordered products to consider. The two local matrix elements are 
given in terms of p\ and p 2 , while the time ordered products are parameterized according to 



-1) 



d\ |T[/C (1) (0)/C 



d 4 x (H(v)\T\lC {1 \0)g {1 



x 



\H(v)} = -27\Tr [m(v)M(v)} 
\H(v)} = -2T 2 Tr { lxTo M(v)s x M(v)} 



d 4 x(iy(u)|T[^( )(0)^W(x)] \H(v)) = -Tr {T al3 M(v){s a , s p }M{v)} 



where 



T 



1 

-a/3- 3 

Using this parametrization one obtains 



Ti '{v a vp - g afi ) + -T£ 'e^pxv^l 7s 



(51) 
(52) 
(53) 

(54) 



m H = M H 



2m, 



■(\i + d H \ 2 ) 



1 



2m c 



[ Pl + 2T X + 2T 3 (1) + d H (p 2 + 2T 2 )] + 0{l/m%) (55) 



where the spins symmetry breaking contribution of the double insertion of the chromomag- 
netic moment operator does not contribute, since we are dealing with the flavor diagonal 
case. 

The parameter pi has been estimated above; this term contribute only about 0.5 MeV to 
the mass of the B meson. The parameter p 2 is more difficult to estimate, but a reasonable 
guess is certainly \p 2 \ ~ motivated by the fact that |Ai| ~ | A2 1 - 

The time ordered products are much harder to estimate. They require in general a model 
describing the dynamics of the light degrees of freedom. We shall not consider this here, 
but it seems reasonable that the time ordered products are of similar magnitude than as the 
local terms. 
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3 Higher Order Corrections to the V c b Determination 



As the main application we consider the higher order corrections to the semileptonic decay 
of a B meson into a D or D* meson. In this case we have to deal with to heavy flavors b and 
c, and the corresponding static operators are denoted b v and c v respectively. These higher 
order terms have been considered already in [|l^] also off the non- recoil point; however at the 
point v = v ' the analysis simplifies drastically compared to the one off the non-recoil point. 

The form factors to be considered are the ones of the vector and the axial vector current, 
defined by 

(fl(p)|&7„c|£>(p')) = yJm B m D h+tvv')2 Vfi + ■■■ (56) 
(B{p)\bjfj,j 5 c\D*{p',e)) = ^m B m D , h A i(vv')(l + vv')e^ H (57) 

where the ellipses denote terms which vanish as v — > v' due to their kinematic prefactors. 
Here b and c are the fields of full QCD and \B(p)) and \D(p')) are the full QCD states. Both 
form factors h + and h A \ are normalized at the non-recoil point v — v' in the heavy quark 
limit such that h + = h A i = 1. In addition to these we also consider the matrix element 

( J B*(p, e )|6 7M c| J D*(p , , e / )) = y /m B .m D .h 1 (vv')(-ee')v, 1 + ■ ■ ■ (58) 

which we shall need to derive normalization conditions. 

Using the 1/toq expansion ([]]) for both operators b and c the contributions to the matrix 
element at the non-recoil point may be classified into three species 

(H b (p)\bTc\H c (p'))\ v=v , = L + T + M + 0(l/m 3 c ) (59) 

where T = 7^,7^75 and Hb and H c are B, B* or D, D* respectively. 

The contribution L are all local terms, i.e. the ones which do not contain any time 
ordered product. They originate from the expansion of the operators ([!]) and read 

L = (H b (v)\b v Tc v \H c (v)} (60) 

-L) {H b {v)\b v T{i^)c v \H c {v)) - Q-L) (H b (v)\b v (z^ 1 ) Tc v \H c (v)) 

^-) 2 (H b (v)\b v T(ivD)(iI? x )c v \D(v)) - (H b (v)\b v (U? L )(wD) Tc v \H c (v)) 



-) (H b (v)\b v (i^ 1 ) Y{i^)c v \H c {v)) + 0(l/m 3 ) 



Am b m c , 

where now \H b (v)) and \H c (v)) are the states in the infinite mass limit. 
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From the discussion of the forward matrix elements it follows that only the last term does 
not vanish. The terms of first order in I/tuq are forward matrix elements of a dimension four 
operator and hence zero, the terms of order l/m b and 1/m 2 vanish after a partial integration, 
which for the forward matrix elements does not yield a surface term. Only the mixed term 
of order l/(m b m c ) yields a contribution, which may be related to Ai and A2. 

The second class of terms are the time ordered products of the current to leading order 
with the terms of order 1/m and 1 / m 2 of the Lagrangian. One obtains 



T 



+ 
+ 



+ 
+ 



-i) 



( 1 



\2m c 
( 1 



d A x(H h (v)\T 



b v Tc v C { c\x)} \H c {v)) 



(61) 



-i) 



\2mi 
( 1 



d 4 x (H b (v)\T [C^(x)b v Tc v \ \H e (v)) 
d 4 x (H b (v)\T[b v Tc v Cf\x)] \H c (v)) 
d 4 x (H b (v)\T [C { b \x)b v Tc v ] \H c (v)) 
d 4 xd 4 y (H b (v)\T [b v Tc v C^(x)£^(y)] \H c (v)) 



\2m c 
( 1 



\2mi 



-if 


( 1 V 


2 


\2mJ 


-*? 


( 1 )' 


2 


\2mJ 




f 1 


-o'l 


\4m b m c 



2 {2^) J xd y {Hb ^ T [Q '^)Q '(y) b ^v\ \h c (v)) 

d 4 xd 4 y (H b (v)\T \c£\x)KTc v C { < f ) {y)\ \H c (v)) 



where here and in the following we suppress the argument of the current b v Tc v which is 
x = 0. 

Finally, there are the mixed contributions M containing a first order term of the expansion 
of the operators ([I]) and a first order term of the Lagrangian 



M 



+ 



2m c 
1 

2m b 
1 



d 4 x(H b (v)\T iKTit^cX^ix)] \H c {v)) 



(62) 



d 4 x(H b (v)\T 



C { b \x)b v (t^) Tc v 



\H c (v)) 



Am c m b 
1 

Am b m c 



d 4 x(H b (v)\T 4 1) (^)^r(z^ i )c„ \H c {v)) 



d 4 x(H b (v)\T 



b v (i^) TcX^i 



\H c (v)) 
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As discussed in section 2.2, these mixed terms all vanish due to the equations of motion. 

In order to proceed further with the time ordered products one has to split the La- 
grangians C h , into its kinetic and magnetic terms. In this way one may analyze the spin 
structure of terms involving products of chromomagnetic moment operators by employing 
the trace formalism, and by using the fact that all products of Dirac matrices may be reduced 
using the algebra of the Pauli matrices, eq. fllOD . The trace formalism gives for the terms of 
order 1/m 



-i) d A x(H b (v)\T 



b v Tc v K,V{x)] \H c {v)) = -xx Tr {m(v)TM(v)} 



(63) 



-i) 



-i) 



d A x(H b (v)\T [b v Tc v gW{x)] \H c (v)) = -xs Tr { lxl5 M{v)Ts x M(v)} 

\H c {v)) = -xx Tr {M(v)TM(v)} 
\H c (v)) = - X 3 Tr {w 5 M(v)s x TM(v)} 



d A x(H b (v)\T 



b v Yc v fC b 



d A x (H b (v)\T b v Yc v Q { 



Here only two parameters xi an d X3 appear since the matrix element has to be symmetric 
under the exchange of b and c and the corresponding exchange of initial and final state. The 
spin structure of the second order terms of the Lagrangian is the same and one may write 
in a similar fashion 



J d 4 x(H b (v)\T[b v Tc v IC^(x)] \H c (v)) = -Sj Tr {M(v)VM(v)} 



(64) 



-i) 



-i) 



d A x (H b {v)\T 
d 4 x (H h (v)\T 



b v Tc v Gf\x)\ \H c {v)) 



-E 3 Tr { lxl5 M(v)Ts x M(v)} 
-Hi Tr [M(v)YM(v)) 
d 4 x(H b (v)\T[b v Tc v gl 2 \x)} \H c (v)) = -S 3 Tr { lx75 M(v)s x TM(v)} 



b v Tc v fC b 



(2), 



\H c (v)) 



Finally, the double insertions of the first order terms are parametrized by 



-i 



d 4 xd 4 y(H b (v)\T [b v Tc v JC^(x)JC^(y)] \H e (v)) = -A Tr [m(v)TM(v)) (65) 
d 4 xd 4 y(H b (v)\T[b v Tc v IC^(x)g^(y)] \H c (v)) = -B Tr { lxl5 M(v)Ts x M(v)} 

d A xd 4 y (H b (v)\T[b v Tc v gi 1 \x)gi 1 \y)] \H c (v)) = - Tr {C a(3 M(v)Ts aS/3 M(v)} 
= -d Tr [m{v)YM{v)) - C 3 Tr [ lxl5 M(v)Ts x M(v)} 
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where we have defined 



1 i 
C a/ 3 = -Ci{-g a p + v a v p ) + -C 2 e Ata/3A w M 7 A 75 



(66) 



A similar expression is obtained for the double insertion of the first order Lagrangian of 
the b quark, involving the same parameters A, B, C\ and C3 due to the exchange symmetry 
b «-> c. 

Finally, the mixed double insertions need another set of parameters 



-1) 



-1) 



-1 



d 4 xd 4 y (H b (v)\T [jC^{x)b v rc v JCP(y)] \H e (v)) = —D Tr {M{v)TM{v)} (67) 
d 4 xd A y{H b {v)\T[K,i 1 \x)b v Tc v gi 1 \y)] \H c (v)) = —E Tr { lxl5 M(v)Ts x M(v)} 
d 4 xd^y{H b {v)\T[gl 1 \x)b v Va u K < i\y)] \H c (v)) = -E Tr { lxl ,M(v)s x TM(v)} 
d*xd*y(H b (v)\T[gi 1 \x)b v rc v gP(y)] \H e (v)) = - Tr {li^M(v)s a Ts p M(v)} 



where R is given in terms of two parameters 



1 i 
R af3 = -Ri(-g af 3 + v a vp) + -i? 2 WA^7 A 75 



(68) 



3.1 The 



and 1 



1 Vector Current at Zero Recoil 



In order to obtain the vector current, i.e. the two form factors h + and hi, we set T = 7^. 
We shall discuss the case of two 0~ states keeping the parameter du = 3 explicit. The form 
factor hi may then be obtained by setting du = — 1 and by replacing the masses mg — > mg. 
and m D — > m D *. 

The matrix element of the local term originating from the expansion of the current is 



(B(v)\b v lli (iI} ± )c v \D{v)) = 2v^M B M D (Xi + d H \ 2 ) 

The traces become trivial and one obtains for (EST) 



(B(p)\hAD(p'))\v. 



( 1 



2JM B M D v„ 1 



' + 



\2m, 



2m b 



\2m c 2m b 
i + A + Ci + d H [E 3 + B + C 3 ]) 



(Xi + d H Xs) 



(69) 



(70) 
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In order to extract the form factor h + from this one has to take into account another 
trivial source of 1/m corrections, which is the normalization of the states. The right hand 
side of Q55| ) is expressed in terms of the physical meson masses m B and m D . In flTU| ) only the 
masses of the static limit appear which differ from the physical masses at the order 1/m 2 . 
To this end, one has to take into account a factor 



IM B M D 1 , ( 1 1 \ 1 
i\ 2m 2 b 



when extracting h + from ([70|) . 

The parameters appearing in flTU| ) are not all independent. The normalization of the flavor 
diagonal current is known in full QCD for both matrix elements 0~ — > 0~ and 1~ — > 1~ 

(5(p)|&V|£(p)> = 2m^ = (B*(p,e)\h,b\B*(p,e)) (72) 

This may be employed to obtain relations between the parameters. Setting mj = m c , ([72]) 
implies the relations 

Xi = X3 = (73) 
2(H 1 + A + Ci + A 1 ) = -(D + i2 1 -Ai) (74) 
2(~ 3 + B + C 3 + A 2 ) = -(2E + i? 2 - A 2 ) 

The first of these equations is Lukes theorem, stating that there are no first order corrections 
at the non-recoil point |TTJ . Taking the relations between the parameters of the second order 
into account, one obtains for the form factor h+ at the non-recoil point 

Ml) = 1 " - \ (D + Ri-^i + 3[25 + R2 - A 2 ]) + 0(l/m 3 ) (75) 

Similarly, by the appropriate replacements one obtains 

hlil) = 1 " (^T " \ {D + Rl ~ Al " [2E + R2 ~ X2]) + °^ m ^ ( 76 ) 

Looking at the definition of the parameters entering the 1/m 2 corrections it turns out 
that to order 1/m 2 the only input need are the two parameters Ai and A 2 from the local 
dimension 5 operators and the time ordered product involving insertions of both and 
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which is given in terms of four parameters. The results for h + and h\ may also be 
written as 



^= i -(^-2s) a K- Ai - 3Aa (77) 

1 r d 4 xd 4 y(B(v)\T\C i b 1 \x)b v c v C^(y)} \D(v))) + 0(l/m 3 ) 



+ -» 



= 1- 



2v/M B M D 
/ 1 1 



+ -i 



\2m c 2m& / 2 
1 



2^M B M D 



-X1 + X2 

d A x d 4 y (B*(v,e)\T (x)b vCv C^ (y) 



(7s; 

I £>>,e)) J +0(l/m 3 ) 



This relation has a simple intuitive interpretation. The contributions from the local dimen- 
sion five operators (Ai and A2) originate from the matching of the field operators of full 
QCD to the ones of the effective theory to order l/m c and l/m& respectively. However, also 
the states receive corrections and the matrix element involving the time ordered product 
corresponds to the corrections to the states to order 1/my, and l/m c respectively. The local 
and the nonlocal contributions are of order 1/ {Am c m^), but due to the normalization of the 
matrix element for the flavor diagonal case all other terms of order 1 / m\ and 1 / m\ have to 
be related to these such that the normalization is preserved in the case m c = m^. This fact 
leads to the prefactor (l/m c — 1/rrib) 2 in front of the correction term of h + (l). 

This result agrees with the one found in In particular, one may see that at the 

non-recoil point the form factors may be expressed in terms of the parameters D\, D3, D4, 
D§ and Dq defined in [I2|. However, the representation in terms of the Pauli matrices reveals 
a relation between the parameters of [I2|, at least at the non-recoil point. In total there 
are six independent parameters Ai, A2, D, E Ri and R2 at v — v' and one may show that 
Ri = S(D 4 + D 5 ) and R 2 = 2(D 5 + D 6 ). 

The 1 / m?Q corrections to h + and h\ depend on the spin symmetry conserving contribution 
X = D + Ri — Ai and on the spin symmetry breaking combination Y = 2E + R2 — A2. If one 
considers in addition the form factor Hax then a third combination of the six parameters Ai, 
A 2 , D, E Ri and R 2 is needed. 

In order to perform a model independent extraction of V c b from the decay B — > Div 
one has to take into account also the second form factor h- of the vector current. However, 
considering only forward matrix elements one cannot say anything about this form factor, 
and one has to consider different velocities along the lines of [12[ . 
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3.2 The — ► 1 Axial Vector Current at Zero Recoil 

For a model independent extraction of V c b the process B — > D*eu is much more interesting 
than B — > /}£z/. The relevant form factor is /i^i as defined in ([57]). To leading order we have 
at the non-recoil point due to spin symmetry 

(B(v)\b v c v \D(v)) = 2^M B M D = -(B(v)\b v (se)c v \D*(v,e)), (79) 

and thus Hai is normalized in the same way as h + and h\. To subleading order Hai will 
differ from both h + and h\ due to spin symmetry breaking. This is, however, calculable in 
terms of the parameters which have been introduced above. 

The local dimension five terms may be expressed in terms of Ai and A2 

(B(v)\b v (iB- L ) llll5 (tIJ ± )c v \D*(v,6)) = 2y/M B M D e fl (~Xi - A 2 ) (80) 

while the contributions from the T-products are evaluated by replacing F — > 7^75 and using 
the representation matrix for the vector meson in the final state. Taking into account the 
contributions form the normalization of the states one obtains for the form factor Hai at 
zero recoil 

Ml) = 1 - (^-) i (£> + i?i - A, + 3[2£ + fl 2 - A 2 ]) (81) 

-{^}\ (D + Rl - Xl - [2E + R ^ H) 

+ (sk) (° + 2E - - R > + + A2 ) + ° (1/mS) 

The structure of this result may be understood from spin symmetry. In the heavy quark 
limit, spin symmetry relates all the form factors h + , h\ and Hai- If one would take into 
account only the corrections of order l/m£, then the spin symmetry of the c quark would 
still be unbroken and one may rotate the D* meson into a D meson. Hence the l/mf 
corrections to Hai have to be the same as the I /ml corrections to h + . Similarly, and more 
importantly, the l/m^ corrections to Hai have to be the same as the corrections to 

h\ since we may now use the spin symmetry of the b quark to rotate the B meson into a 
B* meson. Finally, the mixed insertions break both spin symmetries and thus cannot be 
expressed in terms of hi or h + . 

In total, the three form factors may be reexpressed in terms of three parameters X, Y 
and Z 

( 1 1 \ 2 1 

h + = l-( t:(X + 3Y) (82) 

\m b m c J 2 
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'»= 1 -(^-^)>- y > 2 < 83) 

k - = 1 - ¥ x + 3y > - (i:) 2 V x - y > + s^k (-1* - y + z ) < 84 » 

where 

X = D + R x - Ai, Y = 2E + R 3 - A 2 , and Z = ~D + (85) 

o 

where X corresponds to the spin symmetry conserving interactions, Y to the spin symmetry 
breaking ones, while Z is a mixture of spin symmetry conserving and spin symmetry breaking 
terms. 



3.3 Discussion of the Results and Quantitative Estimates 

Finally we shall discuss the results and try to give a numerical estimate of the corrections. 
In general, this needs input beyond heavy quark effective theory as e.g. a model. A few 
things, however, may be said with two plausible assumptions. 

The form factors h + and h\ are form factors which are in the heavy quark limit related to 
matrix elements of conserved currents, which generate heavy flavor symmetry. The operators 

K + = J d 3 x b v c v , K_ = J d 3 x c v b v , K = J d 3 x [b v b v — c v c v ] (86) 

are generators of the heavy flavor symmetry satisfying [K + , KJ\ = K , and hence one 
derives in the symmetry limit (B(v)\Q+\D(v)) = \JMbMd and a similar relation for two 
vector mesons, implying that h + (l) = hi(l) = 1. 

In the presence of explicit symmetry breaking one splits the Hamiltonian in a symmetry 
conserving piece H and a symmetry breaking term XH S b, such that [Kj , H ] = 0. Since 
now [Kj , H] = X[Kj , Hsb] the generators Kj become time dependent and one has at 
x = 

l= (B{v)\Q + \D(v)) 
2^M B M D 

where X(v) is a hadronic state in which the c quark moves with velocity v. On the left hand 
side we have neglected the matrix element of c v c v between the B meson states. 

Relation ( |§7D is the standard derivation of the Ademollo Gatto theorem [^, and it 
allows two observations. The state X is not in the lowest spin symmetry doublet and hence 
the energy difference Eb — Ex is not of the order A, but of the order 1. Hence the corrections 




(B(v)\Q + \X(v)) 



2^M B M_ 



x 



17) 
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due to symmetry breaking are of second order in the symmetry breaking interaction, which 



is the well known Ademollo Gatto theorem 23 



Secondly, and more importantly for the present discussion, it shows that one expects that 

(B(v)\Q + \D(v)) 



2M B M D 



< 1, 



since the sum on the right hand side of fl8"T|) is positive. This means that h + — 1 < and 
hi - 1 < 0. 

However, it is known that short distance contributions may change this conclusion. For 
instance, the full one loop QCD calculation yields for h + (l) ]25[ 



Ml) = l + -f m6 + mC l*f^)-2) (88) 



7r \rrif, — m c \m, 

which yields a positive contribution to the normalization. This may be traced back to the 
matrix element of c v c v which we have neglected on the left hand side of fl8T|). We shall 
assume in the following that the positive short distance contribution (|38D is compensated by 
the long distance one, for which fl3?D holds. 

From this requirement one obtains two constraints for the parameters 

(D + Ri) - Ai + 3(2E + R 2 ) - 3A 2 > (89) 
(D + R 1 ) - Ax - (2E + R 2 ) + A 2 > 0, (90) 

which are equivalent to (D + Ri) > A x and -(D+Rx-X^/S < (2E + R 2 X 2 ) < (D+R 1 -\ 1 ). 

In fig.|I] we plot the spin symmetry conserving contribution of the time ordered products 
(D + Rx) versus the spin symmetry breaking part (2E + R 2 ). The allowed region is the one 
below the dashed and above the solid line, where hi — 1 < and h + — 1 < 0. 

In order to obtain some numerical estimate we shall assume that (D + R±) < — Ai which 
should be a reasonable order of magnitude for the spin symmetry conserving terms of the 
time ordered products. Thus the parameters for the time ordered product terms should lie 
within the shaded triangular region in fig.|T]. 

We shall estimate the contributions to /iAi(l) by observing that, numerically, the contri- 
butions of the order l/m 2 c are by far the largest. As argued above, spin symmetry enforces 
that these contributions are the same as the ones to the form factor hi. Maximizing the 
form factor h + in the shaded triangle of fig.JT], we have for /^(l) 

4 Ai) <M1)-1<0 (91) 



2m c J \ 3 
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+ 



-0.2 



hi = -0.3 GeV 2 
?i 2 = 0.12GeV 2 




-0.4 



-0.2 



0.2 



0.4 



(D + Rj) {GeV 1 ) 



Figure 1: Allowed region for the spin symmetry conserving 
(D + R±) and spin symmetry breaking terms (2E + R2) of 
the time ordered product. The solid (dashed) line is from the 
constraint h + (l) < 1 (hi(l) < 1), while the vertical dotted 
line is the assumed upper limit for (D + Ri). 



Using Ai ~ —0.3, corresponding to the minimal value assumed here, one obtains corrections 
to hjii — 1 ranging between zero and -5%. This is consistent with the estimate performed in 



12] based on a simple wave function overlap model, once updated values for the parameters 



are used [ [2~6[ 



The present estimate is not based on a model but on the assumption that the spin 
symmetry conserving contributions to the time ordered products satisfy Ai < (D + Ri) < 
— a\\ with a ~ 0{1). The dependence on a is not very strong; for a = 2 one obtains 
Iia\ — 1 > —6% and 1/m^ corrections to h^i exceeding 8% are very unlikely. 



22 



4 Conclusions 



In this paper we have considered forward matrix elements of local operators of higher dimen- 
sion and their time ordered products with terms originating form the heavy mass expansion 
of the Lagrangian. Due to the projection P + = {jj + l)/2 appearing in heavy quark effective 
theory the Dirac algebra simplifies and only two types of matrix elements of local operators 
appear. In addition, the spin structure of the time ordered products of these operators with 
higher order terms form the Lagrangian may be analyzed in a simple way. 

Matrix elements of this type appear in two important applications. Performing a heavy 
mass expansion for inclusive decays along the lines of Bigi et al. ]14[ these matrix elements 



parametrize the non-perturbative input required beyond the leading order in the 1/rriQ 
expansion of total rates as well as for inclusive decay spectra. 

The second application are the form factors for weak transitions at the non-recoil point. 
The symmetries of the heavy quark limit yield the normalization of the weak transitions 
between heavy quarks; this fact may be employed to perform a model independent determi- 
nation of \Vcb\- The recoil corrections to the normalization are given in terms of the forward 
matrix elements considered here. At the non-recoil point the analysis simplifies drastically, 
mainly due to the simplification of the Dirac algebra, as compared to the general analysis. 

As an example we have reconsidered the second order corrections to the semileptonic 
transition B — > D". These corrections have been studied already in |T3 for the general 



case. At the non-recoil point the present analysis agrees with the one performed in ||12j| . 
However, it turns out that at v — v' some of the parameters given in [jn| are in fact not 
independent. 

The second order corrections of to the weak decay form factors are all parametrized in 
terms of five matrix elements 

Ai = ^-(H(v)\Q v (iD) 2 Q v \H(v)), A 2 = ^(H{v)\Q v {iD a ){iD p ){-ia aP )Q v \H{v)) 
and the matrix elements of double insertions of the first order correction to Lagrangian 
Ml 1 



-%. 



2y/M B M D 
1 

2y/M B M D 
1 



d A xd 4 y (B(v)\T [K$\x)b v c v 1CM{yj\ \D(v)) 
d 4 xd 4 y (B(v)\T [gPixfocvKPiy)] \D(v)) 
d 4 xd 4 y (B(v)\T \Q£\x)b v c v g^{y)\ \D(v)) 



2^M B M D 

All other matrix elements of time ordered products are related to these by the normalization 
condition for the vector current in the full theory. 
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All these matrix elements are non-perturbative. In principle they may be measured on 



the lattice and first results have been reported p7| . However, in the meantime one has to 
rely on some model to estimate their size. In the present paper we have used a reasonable 
guess for the spin symmetry conserving contributions of the time ordered products to get 
some upper limit for the I/tuq corrections to h^i at the non-recoil point. The main result 
of this analysis is that the corrections to /iAi(l) are small, ranging between —5% and zero. 



Including also the leading and subleading QCD radiative corrections |10[ to the normal- 
ization of Ua\ one concludes that 

= 0.96 ±0.03, (92) 
and thus the corrections at zero recoil are small. 
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